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Propagation of Very-Low-Frequency Pulses to 

Great Distances 1 

James R. Wait 

A theoretical study is presented for the propagation of electromagnetic pulses at very 
low frequencies to large distances. The space between the earth and the ionosphere is 
represented as a wave guide with sharply bounded and concentric spherical boundaries. 
The concept of phase and group velocity and its application to the present problem is dis- 
cussed in some detail. The influence of the propagation medium on the shape of the envelope 
of a quasi-monochromatic pulse is also considered. Using an alternative approach, the 
response of an impulsive source is also calculated and is shown to be a damped oscillatory 
function of time with a quasi-half-period varying in a predictable maimer with distance of 
travel in agreement with the observations of Norinder and Hepburn. 

1. Introduction 

In recent years there has been a renewed interest in the propagation of VLF (very-low- 
frequency) radio waves. The wave-guide model has been remarkably successful for the pre- 
diction of the field strength as a function of distance and frequency. In most investigations, 
both experimental and theoretical, the source varies essentially in a sinusoidal manner with 
time. 2 It is of great interest, however, to have a knowledge of the transient characteristics 
of a pulse and the way it is influenced by the propagation medium. For example, if the re- 
ceived field strength of a lightning stroke is recorded as a function of time, one may ask how 
this waveform is related to the source current, distance of travel, electrical constants of the 
ground, and height and electrical properties of the ionosphere. It is the purpose of this paper 
to attempt an answer to this question on the basis of a theoretical study. 

A general analysis is presented in section 4, applying and extending the classical concepts 3>4 
such as phase and group velocity for propagation in dispersive media. Such an approach is 
particulaily suitable for a quasi- monochromatic source wherein the spectral components are 
centered in a narrow band about some central or reference frequency. In section 5, an alter- 
native method is described that is applicable to broad-band sources which contain many spectral 
components. Extensive numerical results are presented for the parameters characterizing the 
shape of the radiated pulses. Finally, some reference is made to published experimental data 
concerning the waveforms of radio atmospherics. 

2. Quasi-Monochromatic Pulse 

In the present study, the earth is represented by a homogeneous sphere of radius a, and the 
ionosphere is idealized as a sharply bounded concentric reflecting layer at height h. The source 
of the field is always assumed to be equivalent to a vertical (radially oriented) electric dipole. 
At large ranges, the field can then be represented as a sum of waveguide modes. 

The vertical electric field of any one mode at a range d can be written in the form (see sec- 
tion 7.1) 

Jbcd .2ird 

e(co)=A n e hUn e ' x a V"S (1) 

where A n is a slowly varying function of frequency co, u n is a measure of the attenuation of the 



1 The results in this paper were i eported at the commission 4 and 6 sessions of the International Scientific Radio Union General Assembly held 
in Boulder, Colo., August 22 to September 5, 1957. 

2 See VLF issue of Proc. Inst. Radio Engrs. 45, June 1957, for papers by Budden, Pierce, Wait, and Watt. 

3 J. A. Stratton, Electromagnetic theory, p. 292 (McGraw-Hill Publishing Co., Inc., New York, N. Y., 1941). 

4 IA. L. Al'pert, V. L. Ginzburg, and E. L. Feinberg, Radiowave propagation, pt. II. p. 364 (Moscow, 1953). 
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mode per unit distance, and s n is a dimensionless phase factor. The phase velocity is then 

given by 

v p =c/s n , (2) 



where c is the velocity of light. On the other hand, the group velocity of a mode is given by 

c 



g ~~d 
which can be written in terms of H(=h/\) as follows: 



(3) 



o 4-H- n 



(4) 



dH 



The group velocity, so defined, is a measure of the velocity of the envelope of the quasi-mono- 
chromatic pulse or group. For example, if the propagation medium were nondispersive, the 
signal might be represented by the real part of E (t) =A(t)e iu * 1 , where co is the carrier frequency, 
and A(t) is the shape of the envelope as a function of time. In the simple case of a broken 
sinusoid, 

A(t)=l for -f<*<f» 

A(t)=0 for t>T/2 and t<-~ 
The group velocity concept is valid only when the frequency spectrum, 

G(o>)= f + ° ' E {t)e- iat dt, (5) 

%) — oo 

differs substantially from zero only in a small frequency band near the signal carrier-frequency 
w . In the present case 

G( u )=2 8in( r^ T ' 2 - (6) 

(co— OJ ) 

The spectral width is given by Aco^l/(2T) and this must be small compared to co . In other 
words, the period of the envelope must be long compared to the period of the carrier frequency. 
Now when the pulse propagates through the medium by a distance d, the field E (t) of a 
mode becomes transformed to E(t) whence 



*®-sZ 



AWe^'-tWe-^Gicdda, (7) 

an j _ a, 

where 

and 

, N 2ird 

The coefficients a(co) and </>(o>) are phase and attenuation factors of a mode as a function 
of frequency for the section of path of length d. 

Utilizing well-known properties of Fourier integrals, 1 2 the spectral representations of the 
signals E (t) and E(t) can be written 

#<>(*) =;r A(t')exp{i(ca —u)t'+a>t}d<adt', (8) 

ATTJ- oo J-oo 
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and 



£(<)=J-f "f C °^(Oexp{i(« -co)^+^-^(<o)}G'(a,)e-«<»V/corf^. 

^"%J — 00 %) — oo 

The phase function $(co) is now written in a Taylor expansion, 



(9) 



O 2 

*(«) = ^(co ) + o</>' (o>„) + T <t>" («»)+■.. , (io) 

where Q=o)—co . Neglecting terms containing O 3 , etc., it follows that 

~r^ J + °° A(t')exp SiQlt-tr-i'^-iQ-^&Xdl'dQ, (11) 

where the slowly varying factors have been taken outside the integrand. Introducing a new 
variable z defined by 

2 =V?W»[«+^±^} (.2) 

in place of $2 it follows that 

E(t) =~^ exp -J i Lc () f - 4>(u ) 1 V G(u ) e— W> X 



since 



exp — i-z 2 \dz=l—i. (14) 

After a further change of variable, via 



4>»M ~^ U > (15) 

we finally arrive at 

E(t)=~ exp{i[^t-^ Q )]}G^ Q )e-^)xj^ A[t-<$>' {u,)-^' l^)u] exp (i| u^j du. (16) 

In the above, it has been assumed that 0"(coo)>O; if <l>" (co ) <^0 , then by changing the sign in 
the substitution of variables one would have -yVl^'^o)! i R pl& c e of ^1 T( f)" (a> ) . 

Now if // (w o )=O and noting that 

f exp(i^vAdu==l + i, (17) 

it follows that 

^(^=4^-0' (woJle-^MGCcooJe^-*^)]. (18) 

To this approximation, the pulse has not changed shape, although the pulse amplitude has been 
modified, the phase has changed by an amount <£(co ), and the puls e as a whole is delayed by the 
"group time" </>'(W). The effect of a nonzero value of V W" (w ) is to modify the pulse shape 
(see footnote 4). To evaluate this phenomenon, the envelope must be specified and the inte- 
gration indicated in eq>(16) must be carried out. For example, in the case of the broken 
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sinusoid extending from — T/2 to T/2, the integral in eq (16) can readily be transformed to a 
Fresnel integral, 



E(t)=-j- exp {i[a> t— cj>(coo)]} I exp \i^u 2 \du, 



(19) 



where 



u x = 



vW'W 



u 2 = 



T-e 



V**"(«b) 



and 0=(T/2)+£— 0'(w o ). is the time measured from the instant — T/2 + <t>' («o) • In the case 
of negligible dispersion of the pulse shape, 0"(o> o ) must be sufficiently small to enable u x and u 2 
to be replaced by — oo and +°°, and then 



E(t)o*exv{i[a> t-<t>(uo))} for O<0<T 

^0 for 0<O and 0>T. 



(20) 



To illustrate the influence of dispersion, it will be assumed that the signal duration T is large 
compared to V 't4>" (co ) . Then the form of the leading edge of the signal is determined by 



™i^ir e 4*HH 



(21) 



The envelope \E(t)\ is plotted versus the time parameter fl/vW" (co ) in figure 1, a. It is now 
convenient to define a build-up parameter t b , which is the time measured from 0=0 for the 
envelope \E(t) | to approach within 5 percent of unity. From inspection of figure 1 it is seen that 



4^4-v / 7r0' / (coo) sec. 
For numerical presentation it is convenient to introduce a dimensionless build-up parameter 



T b defined bv 



where 



4<s/dh/2 



T» 



(22) 



T b = 



ds n -^d 2 ^ 



The build-up parameter T b and the build-up time t b are related by eq (22). The chart in 
figure 1, b, is to facilitate the conversion for a typical height of 70 km, and for various ranges 
d in kilometers. 
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Figure la. 
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3. Numerical Presentation of Waveguide Parameters 

The important parameters in specifying the characteristics of pulse propagation in the 
earth-ionosphere waveguide are attenuation rate, phase velocity, group velocity, and the 
build-up parameter. In general, these are functions of the dimensions of the guide and the 
electrical properties of the bounding walls. Those four quantities (u n , v p /c, v g /c, and T b ) are 
plotted in figures 2 and 3 for mode number 1 (i. e., the dominant mode). The abscissa is 
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H(=h/X) } which is the height of the ionosphere in wavelengths and may be regarded as a 
frequency parameter. 

The dimensionless quantity A is a ground conductivity parameter and is defined by A= 
GjH, where G= e o)/a g in terms of the conductivity cr g of the ground. Alternately, one may write 
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The dimensional quantity B is the corresponding parameter for the (assumed) sharply bounded 
ionosphere. It is defined by B=L/H, where Z— o>/a> r , 

A 2/ (plasma frequency) 2 

and o) r = coov= — yp—. 1-= 

colhsional frequency 

In terms of the effective ionospheric conductivity, a {j 5 B is given by 

p _ 1 
60<r t h 

The curves in figure 2, a to d, inclusive, are for ^4=0 or o- g ^°° . This assumption is valid 
for propagation over sea water for all values of H shown. The curves in figure 3, a to d, in- 
clusive, are for yl = 10~ 4 , which corresponds to moderately conducting ground [2.2 millimhos/m 
for h = 70 km]. 

4. Analysis for Impulse Response 

The foregoing analysis is applicable to the description of wave packets wherein the spectral 
components are essentially contained in a relatively narrow band about the carrier frequency. 
In the case of a lightning discharge, the pulse contains many spectral components, and as a 
matter of illustration, if the current was in the form of a Dirac delta function, all frequencies 
would be contained in a uniform distribution. In this case and for other similar instances, 
it is desirable to use an alternative approach. In the following, it is assumed that the source 
is initiated at t=0 and d=0. 

The most convenient description of the source is the (transient) waveform of the radiation 
component of the vertical field on a perfectly conducting flat ground plane at distance d. This 
field is denoted e (t) as a function of t and its Fourier transform is 



fi 



/» 00 

^ Jo ~ e « (t)e ~ iUtdf - (23) 

If the source is a dipole carrying a current j(t) with a height h(t), both functions of t, then the 
radiation field on the (hypothetical) flat ground plane is 

« oit) =2^\ji ljmm ~] l : l . u{t ' } ' (24) 

where 

u(t') = l for *'>0l 

yt'=t-to,to=d/c. 
t'<o] 

The field at range d of the same dipole source when located on a homogeneous spherical earth 
of radius a with concentric ionospheric reflecting layer at height h is denoted e(t), and its 
transform is 



=0 for V 



E( 



?(w)= r e(f)e- i(at dt. (25) 



The relation between K(co) and Ao(co) has been previously derived (see footnote 2), and the 
final result is simply quoted here. 

EM 4-^rJ ^ #§ ± si- fcT?" .-<•>, (26) 

|_sin d/aj h (Ka) /2 h=o 

where u n and s n are the attenuation and phase constants of the waveguide mode of order n, as 
defined in section 7.1, 8 ^%, and 8 n ^l iorn^O, and F((a) = (u/c)s n d—ut. The waveform of the 
electric field e(t) is now obtained from the inverse Fourier transform 

e(t)=?<r r* ' E^e^'da. (27) 

Z7T J-oo 



5 The concept of an effective isotropic ionospheric conductivity is valid for co«>^10 7 and highly oblique (low-order) modes. 
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Unfortunately, the integration, except in certain limiting cases, cannot be carried out in closed 

form because u n and s n are not simple analytical functions of co. In view of the fact, however, 

that extensive numerical data is already available for u n , s n , and related functions, it seems 

desirable to evaluate the integral by an approximate saddle-point procedure. The integrals 

to cope with are of the form 

*+• 
J= 



= f + " G(^e- iF ^dco. (28) 

t/ — oo 



G(a>) can be assumed to be slowly varying compared to F(w). The saddle point of F(ui) is at 
co=cjo s , which is a solution of 

do) 

Expanding F(u>) in a Taylor series about co s leads to F(a))=F(co s )+[(ca— u s ) 2 /2]F" (w 5 )+ terms 
containing (co— oj s ) 3 , etc., where 



'-w-PB 4 ! 



Retaining only the first two terms in the expansion for F(u>), the integration can be carried 
out to yield 

J "kiF^T*to«--^ (29a) 

where the asterisk denotes a complex conjugate. When the current is in the form of a unit 
impulse or Dirac delta function, the radiated field 'e (t)/is a doublet impulse function. That is, 

e Q (t)=E 1 8'(t) (29b) 

where 

5 (t)=Lim — 

A->0 A 

Therefore, the Fourier transform is 

'E (o))=E 1 i<a. (29c) 

In this important case 

Therefore, the transient response for the impulse current source is given by 

e( ^4^T (2,r K0s 5A( ^ (3i) 

e n (t)= (g,)»[«.K)] 3/, C y t ff exp[=ftao>.)], (32) 



where 
because 



F'\co s )=2 ( ^+ ( L^^[T b (H s )}\ (33) 



The quantities H s , s n (oo s ), u n (u s ), and T b (H s ) are the values of H, s n , u n , and T b evaluated 
at the saddle point a)=co s . 

It can be readily seen that the equation for determining the saddle point is equivalent 
to finding the value of H (or co) that satisfies 

l_v g (H)_d 



(34) 
s n c ct 

Because v g (H)/c is a smooth monotonic function of H, the saddle-point value H s (or o) s ) is readily 
determined by using a simple graphical procedure for any specified value of t. 
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5. Application of the Impulse Responses 

The response of the waveguide to an impulsive current source is shown in figure 4 for various 
ranges. The quantity plotted is e\(t) given by eq (32) with n=L This is the waveform of 
the dominant mode for the typical daytime conditions (k=7() km and L/H=0A) and prop- 
agation over sea water (A^O). It is interesting to note that the pulses have the appearances 
of damped sinusoids. The oscillating nature of the curves is not due to the source which is 
impulsive, but rather is a result of the modal characteristics of propagation medium. In 
general, the initial quasi-period (temporal length of the first half-cycle) is becoming progres- 
sively shortened with increasing range, whereas the oscillatory nature of the pulse is becoming 
enhanced. In each case the frequency of the latter part of the waveform seems to be approach- 
ing about 9 kc (i. e., a half-period of about 56 jusec). 

It is desirable to repeat this calculation for a source current of finite duration. For 
example, if the current dipole moment is proportional to ae~ at then the Fourier transform 
of the primary field e (t) is 

E ( U )=E 1 -^-, (35) 

which replaces eq (29c). The transient response, denoted by e ti (a,t), is thus given by 

a cos [F(o) s )-\-\//(co s )] 



e n (a,t) = (H,Ms n MY 



3/2 . 



(«*+«*)* 



T,,(a s ) 



expl — ^— «,(«,) J> 



0*6) 



where ^(«,)=arctan («.,/«). As a tends to °°, e n (a,t), of course, approaches c N (t), which is the 
impulse response. 6 

The response ei(a,t) for the dominant mode at d=3,000 km is indicated in figure 5 for the 
same daytime conditions as the curves in figure 4. Various values of the time constant 1/a 



i Note that d(t)= -H™ <*«-«' for t>(). 
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are indicated on the curves. The impulse response corresponds to the case 1/a— 0. The 
general nature of the curves is very similar. There appears to be a shift in the phase of the 
cycles but the quasi-half-periods are essentially unchanged. 

A more realistic waveform for the current dipole moment S(t) of a lightning stroke is a 
pulse which rises up smoothly from zero to a peak value in about 10 or 20 jusec and then decays 
to zero somewhat more slowly. Examples of such pulses are shown in figure 6, where S(t) 
is of the form 



S(t) = c 



(37) 



with a t y>a. The response e x (t) to this composite exponential source is then simply obtained 
by superposition of the exponential responses. For example. 



ei(f)=ei(a,f) e x {a u t). 



(38) 



Employing the three particular forms of source current S(f), the response ei(t) for the dominant 
mode is shown plotted in figure 7 for typical daytime conditions over sea water. The curves 
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again are very similar in appearance to the previous response curves for impulse and single 
exponential sources. The effect of finite rise time of the source current is to reduce the ampli- 
tudes of the peaks in the initial part of the waveform. 

It is rather important to note that the quasi-half-periods are not appreciably influenced by 
the nature of the source pulse. It would therefore seem justified to compare the variation of 
the quasi-periods with range for the calculated impulse responses and the experimental wave- 
forms of Hepburn. 7 Such a comparison is shown in figure 8. The solid curves are the calculated 
quasi-half-periods in microseconds for daytime conditions as a function of range. Only the 
first five half-cycles are shown. The corresponding data from Hepburn's paper is shown by 
the encircled numbers which indicate the order of the half-cycle. The agreement is quite 
reasonable. 



7 F. Hepburn, Wave-guide interpretation of atmospheric waveforms, J. Atmospheric and Terres. Phys. 10, 121 (1957). 
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Figure 8. Quasi-half-periods of dominant mode 
(n = l t A^O, h = 70 km, i/H-O.l), 
©From Hepburn's measured waveforms. 

The calculated variation of the quasi-half-period with range and time are also in good quali- 
tative agreement with Norinder's 8 observed waveforms who first pointed out the characteristic 
transformation of the quasi-half-periods. 

The transient waveforms presented in figures 4 to 7 are the dominant (n=l) mode only. 
The same procedure can be used to calculate the higher order mode responses. At ranges 
exceeding about 2,000 km, these do not contribute appreciably to the total field for the range 
of times shown on the curves. 

The transient response of zero-order mode can also be important for longer times. It 
corresponds to what has been called the "slow tail" in atmospheric waveforms. Again for the 
range of time shown in figures 4 to 7 it would not be significant. For the sake of completeness, 
a short analysis of the transient response of the zero-order mode is presented in section 7..'}. 



6. Conclusion 

The analyses and results presented in this paper should be useful in the interpretation 
of experimental results for the propagation of pulses to large distances over the surface of the 
earth. It appears to be desirable to separate the transient analyses into two parts depending 
on whether the source is a quasi-monochromatic pulse containing a narrow band of frequencies 
or an impulsive type having many spectral components over a wide frequency range. The 
parameters describing the propagation of these two classifications of pulses have been presented 
in graphical form. 

7. Appendixes 

7.1. Waveguide Modes 9 

The earth is taken to be a homogeneous sphere of conductivity a and dielectric constant e. 
The lower edge of the assumed homogeneous ionosphere is taken to be at a height h. The 
collisional frequency is v and the plasma frequency is co . 

Assuming that the source is a vertical dipole radiating P kilowatts, the vertical electric 
field in millivolts per meter at the great circle distance d in kilometers, for a time factor exp 
(iwt), is given by the mode sum, 

E=E W : 

8 H. Norinder, The waveform of the electric field in atmospherics, Arkiv for Geophysik 2, 161 (June 1954). 
8 J. R. Wait, On the mode theory of VLF ionospheric propagation, Geofis. pura e Appl. 37, 103 (1957). 
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where 

^ r (d/a) YVWL e *[¥-i]± ss ./v-«-«w (1 1) 

^=Lsin (d/a) J (h/\) e ^o " n ' ( } 

„ 300 VP _r sin(47rC w feA) "|" 1 « I for ^=0 

£-0- d , d B -^i-|- ( iTCnh/X) J ^l for ti^O, 

where (S' n =(l — C£)'^, and C„ is a solution of 

R g (C)R t (C)exp(-4:TiCh/\) = exp(— i2ira). (1.2) 

R S (C) is the Fresnel reflection coefficient at the ground for a plane wave whose (complex) 
angle of incidence is the arc cosine of C, and Rt(C) is the corresponding Fresnel reflection 
coefficient at the lower edge of the ionosphere. 

It follows from the previous analysis (see footnote 9) that 

p ( n, (KO-i)C-[(K-\)Q*-iQ+C*Q*\» n „. 

where K=e/e , and 6=(t^o)/a, and 

where 10 

N 2 ^l-i/L=--l-i(ic r /co), a> T = a>l/v. (1.5) 

The complex values of S n (n=0,l,2, . . .) satisfying eq (2) have been obtained from an 
automatic computor by using a program devised by H. H. Howe. When L is small compared 
with unity, it was shown previously (see footnote 9) that 

° n - L \2hJ J 4irA 6 [ ^ } T 1 /XnVY 2 



nm' 



wiie re 

i =l, e n =2(n5*0). (1.7) 

The above formula for S n is particularly suitable for very low frequencies, in which case 
only the zero mode is of any consequence. 

The exponential term inside the summation of eq (1.1) determines, in the main, the 
propagation characteristics of the modes. It can be rewritten as follows: 

exp[— i2irS n (d/\)] = ex-p — j- u n exp l T s n I (1.8) 

where u n = — ImS n (h/\) J and s n =ReS n . 

7.2. Transient Response of Ideal Waveguides 

Most of the complications in the waveguide mode theory of VLF ionospheric propagation 
are due to the finite conductivity of the bounding walls. The characteristics of the modes are 
obtained only from a numerical iteration procedure. When the losses in the walls are negligible, 
it is possible to derive somewhat simpler forms for the transient response and associated 
parameters. 

Assuming that both walls of the waveguide are of perfect conductivity, the phase velocity 
is given by 



H<H&)T 



!o These expressions for the effective refractive index which tacitly neglect the earth's magnetic field arc valid for VLF and small values of C. 
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and the group velocity by 
Furthermore, 



and, consequently, the saddle point is obtained from 



bH 



0, (2.4) 



which is satisfied by 

where 
Consequently, 



#=# s =|f, (2.5) 



X=[(ct) 2 -d 2 ]'<K (2.0) 



F(HJ—%£, (2.7) 



and 



"s n =s n (H s )=d/cL (2.9) 

Inserting these values into the saddle-point formula given by eq (32) leads to 

for n=l, 2, 3. . . . 

The preceding formula is the waveform for the vertical electric field at range d on an ideally 
perfectly conducting earth with a perfectly conducting ionosphere at height h. It is valid for 
large distances such that d>>/^ and X, 

If the upper boundary was a perfect magnetic conductor (phase shift of 180° on reflection 
rather than 0°), the preceding formulas are modified by replacing n by n — J. 

7.3. Transient Response of the Zero-Order Mode 

The saddle-point method is valid only if the function exp[— iF(<a)] is rapidly varying 
compared to other factors in the integrand of the inverse Fourier integral. For mode zero this 
is not so, but in this case it is known (see footnote 9) that 

where 

a g ^-Jcr e fjf2h and a, ;^V a t fi • 2h. 

Therefore, for the impulse current source, 
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The transient response e (t) of the zero-order mode is thus given by 

e{t)=^ + y(o>)e^"(L>, 
where t' — t—t . This can be evaluated to yield 



where 



and 



with 



d' 






-P/w 



(3.3) 



(3.4) 



(3.5) 



and 



V " V°"i V°* 



The transient response is thus proportional to the characteristic function P (P/t'). In view of 
the combination in which the ionospheric conductivity vi and the ground conductivity a h 
occur, and because (Tg> >(*■*, it follows that a^o^. The transient response of the zero-order 
mode does then not depend to any extent on <r g . The function P(($/f) is plotted in figure 9 as 
a function of T, where T=t'/p. The multiple scale shown at the bottom of the figure is to 
facilitate the conversion of the parameter T to actual time in microseconds. 

The transient response of the zero-order mode has been studied extensively by Schumann. 11 

u W. 0. Schumann, Uber die Oberfelder bei Ausbreitung langer elektrischer Wellen urn die Erde und die Signale des Blitzes, Nuovo cimento: 
9, 1 (December 1952). 
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Figure 9. Zero-mode transient response for impulse source. 
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7.4. Justification of the Saddle-Point Method 

It is worth while to examine the transient calculation from a slightly different viewpoint 
in order to give further justification to the use of the saddle-point method of integration. When 
the earth and the ionosphere are represented by perfectly conducting planes with constant 
separation A, the magnetic field on the earth of a vertical electric dipole on the earth is a sum of 
modes of the type 

H„(«)=^, (4.1) 

where T n is a ^-directed Hertz vector given by 

2iir 2 



*«<Hf[kd(l-^* (4.2) 



« n — 7 £ 

apart from a constant factor. H { 1 } is the Hankel function of the second kind for an impulsive 
current source (i. e., current moment proportional to d(t)), the magnetic field of a mode is then 
given by 

This integral can he evaluated 12 » 13 to yield 

, . v x d c r irnX 2 irnX , ct+X"] ,. AS 

K{t)= iMx L cos ir" sm ~tr log - ~ir} (4 - 4) 

where X=[(ci) 2 — d 2 ] X/ K When X << ct or d and d >> h, the term containing the logarithm 
is negligible, and consequently 

Kit) -^ |cos (-g-) +- r sin (^j-j J- (4-5) 

This latter form could have been obtained directly from eq (4.3), using the saddle-point pro- 
cedure adopted in the main body of the paper. 

7.5. Symbols 
co = (angular) frequency, 

d— distance from source to observer measured along the surface of the earth, 
h = separation between the earth and the lower edge of the ionosphere, 
\= wavelength in free space, 

u n = attenuation factor of the waveguide mode of order n, 
s n = phase factor for a mode of order n, 
c= velocity of light in free space, 
v p — phase velocity of a mode, 
v g = group velocity of a mode, 
H=h,!\ 
E (t) =form of undistorted signal, 
A(t) = envelope of undistorted signal, 

co = carrier frequency of undistorted signal, 
T— duration or rectangular envelope, 
G (co) = Fourier or frequency spectrum of E (t), 
E(t)= form of signal after propagation through dispersive medium, 



1? N. W. McLachlau and P. Humbert, Memorial des sciences mathematiques fascicule 100, 34. Paris (1941). 
1? K G. Budden, The propagation of a radio atmospheric, 42, 1 (January 1951). 



201 



, v 2TT(i 

4>{co) =8 n — — and 

A 

a(w) =u n—r- are phase and attenuation factors of a mode as a function of frequency for a range d, 

/' = an integration variable in eq (8) to (15), 
12 = co — coo, an expansion parameter, 
</>' (co ) =[dc&(co)/dco]co=a> , 

u-i and w 2 =limits of integration in eq (19), 

T 
= the time measured from the instant — ~+</>'(co ), 

/ft = time from 0=0 to the point where the envelope approaches within 5 percent of unity, 

T b = — t b , a dimensionless build-up parameter, 

eo = 8.854 X 10~ 12 , permittivity of free space, 
<r g = conductivity of the ground, 
<7 = €oco/o-£, 
A = G/H=--1/Q0a g h, 
co = (angular) plasma frequency, 
;>= (angular) collisional frequency, 

o-i = effective ionospheric conductivity, 

L = co/co r , 

B=L/H=1/6(W&, 
So (0= the radiated electric field for an ideal flat perfectly conducting ground plane, 
7io(w) =the frequency spectrum of e (t), 
j (0= instantaneous average current of the source dipole, 
h .(/) = instantaneous height of the source dipole, 
5(0= unit impulse or Dirac function, 

M=4xX 10~~ 7 , permeability of free-space, 

h—d/c, 

t =t — to, 

e(t)= resultant field of (transient) dipole sources, 
/?(«)= the frequency spectrum of e(t), 
a = radius of the earth, 
5„^1 for n^O, 

^i for n = (see eq (1) for more accurate definition), 

F(a))=— s n d — cot, 

co s =the saddle point of F(co), 
Ei = constant of proportionality in eq (9b) , 
5' (t)= doublet impulse function, 
T h (H s ) = T h , value of T b at the saddle point H=H s = o3 s h/2irc 7 

C n = cosine of the (complex) angle of incidence of the mode of order n at the bounding walls in the 

waveguide, 
S n =(l -€*)", 
e = dielectric constant of the ground, 
K=e/eo, 

N*=l—i/L=l—iurla, 
7 n = 2 for n t£ 0, 

€0=1, 

X = [(ct)*-d>]x, 



"-am"*®"} 



[J N (a,) = magnetic field of the nth mode of a dipole source in a parallel plate waveguide, 
7r n = the corresponding ^-directed Hertz vector. 
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8. Addendum 

Because of certain normalizations, the distance d between transmitter and receiver is 
measured in terms of the great-circle distance at height h/2. The phase and group velocities 
shown in figures 2 and 3 should then be multiplied by (l—h/2a) ^0.995, if referred to the great- 
circle distance on the earth's surface. In other places in the paper, this correction is negligible 
and can be ignored. 

It should be noted that the correction factor (l—h/2a) can vary because of changes in 
height of the ionosphere and changes in the effective earth radius a. For this reason the veloci- 
ties were presented in figures 2 and 3 in uncorrected form. 

There is a further modification to the attenuation factors when a more recondite mode 
equation is employed. 14 It appears that the attenuation factors are increased by as much as 
50 percent for jFT>4 when ?i=l and ^4=0 but less for other values of n and A. This will not 
change the appearance of the waveform shapes for the range of time indicated in figures 5 to 9. 



The numerical results in section 3 were obtained by H. H. Howe and A. M. Conda, and 
those in section 5 were obtained by W. E. Mientka. Many helpful suggestions were received 
from K. G. Budden, A. G. Jean, E. A. Lewis, K. A. Norton, and W. L. Taylor in the course 
of this work. 



" This equation is 

/,\(0/ ( \(C>Xl> [-47r/r7/(l+^y^]=OX,>(-/L>7r,/), 

which reduces to eq (1.2) for [i?«C] 2 »A/a. 



Boulder, Colo., December 18, 1957. 
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